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Abstract 

We present the results of experimental study of nodal domains of wave functions (electric field 
distributions) lying in the regime of Shnirelman ergodicity in the chaotic half-circular microwave 
rough billiard. Nodal domains are regions where a wave function has a definite sign. The wave 
functions ^'at of the rough billiard were measured up to the level number = 435. In this way the 
dependence of the number of nodal domains H^r on the level number N was found. We show that 
in the limit N ^ oo a least squares fit of the experimental data reveals the asymptotic number of 
nodal domains Hjv/A^ ~ 0.058 ± 0.006 that is close to the theoretical prediction Hat/A ~ 0.062. 
We also found that the distributions of the areas s of nodal domains and their perimeters / have 
power behaviors Ug oc s~'^ and n; oc l~'^' , where scaling exponents are equal to r = 1.99 it 0.14 
and t' = 2.13 it 0.23, respectively. These results are in a good agreement with the predictions 
of percolation theory. Finally, we demonstrate that for higher level numbers N ~ 220 — 435 the 
signed area distribution oscillates around the theoretical limit ~ 0.0386A^^. 
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In recent papers Blum et al. flj and Bogomolny and Schmit have considered the 
distribution of the nodal domains of real wave functions in 2D quantum systems 

(billiards). The condition ?/) = determines a set of nodal lines which separate^egions 
(nodal domains) where a wave function \E'(a;,?/) has opposite signs. Blum et al. [1] have 
shown that the distributions of the number of nodal domains can be used to distinguish 
between systems with integrable and chaotic underlying classical dynamics. In this way 
they provided a new criterion of quantum chaos, which is not directly related to spectral 
statistics. Bogomolny and Schmit have shown that the distribution of nodal domains 
for quantum wave functions of chaotic systems is universal. In order to prove it they have 
proposed a very fruitful, percolationlike, model for description of properties of the nodal 
domains of generic chaotic system. In particular, the model predicts that the distribution of 
the areas s of nodal domains should have power behavior oc where r = 187/91 j^. 

In this paper we present the first experimental investigation of nodal domains of wave 
functions of the chaotic microwave rough billiard. We tested experimentally some of impor- 
tant findings of papers by Blum et al. [1] and Bogomolny and Schmit such as the signed 
area distribution or the dependence of the number of nodal domains Kjv on the level 
number A^. Additionally, we checked the power dependence of nodal domain perimeters I, 
ni oc Z""^', where according to percolation theory the scaling exponent r' = 15/7 [3], which 
was not considered in the above papers. 

In the experiment we used the thin (height h = 8 mm) aluminium cavity in the shape of 
a rough half-circle (Fig. [1]). The microwave cavity simulates the rough quantum billiard due 
to the equivalence between the Schrodinger equation and the Helmholtz equation 4, la]. This 
equivalence remains valid for frequencies less than the cut-off frequency = c/2h ~ 18.7 
GHz, where c is the speed of light. The cavity sidewalls are made of 2 segments. The rough 
segment 1 is described by the radius function R(9) = Rq + Y.m=2 amsin(m6' + (f)m), where 
the mean radius -Ro=20.0 cm, M = 20, and 0m, are uniformly distributed on [0.269,0.297] 
cm and [0,27r]jrespectively, and < < vr. It is worth noting that following our earlier 
experience js, 7] we decided to use a rough half-circular cavity instead of a rough circular 
cavity because in this way we avoided nearly degenerate low-level eigenvalues, which could 
not be possible distinguished in the measurements. As we will see below, a half-circular 
geometry of the cavity was also very suitable in the accurate measurements of the electric 
field distributions inside the bilhard. 
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FIG. 1: Sketch of the chaotic half-circular microwave rough billiard in the xy plane. Dimensions 
are given in cm. The cavity sidewalls are marked by 1 and 2 (see text). Squared wave functions 
l^^NiRc, 0)1'^ were evaluated on a half-circle of fixed radius Rc = 17.5 cm. Billiard's rough boundary 
r is marked with the bold line. 



The surface roughness of a billiard is characterized by the function k{9) = {dR/ dd) / Rq. 
Thus for our billiard we have the angle average k = ((^^(^))5i)^''^ — 0.488. In such a 
billiard the dynamics is diffusive in orbital momentum due to collisions with the rough 
boundary because k is much above the chaos border kc = M~^/^ = 0.00056 ^. The 
roughness parameter k determines also other properties of the billiard 9[. The eigenstates 
are localized for the level number N < = 1/128A;^. Because of a large value of the 
roughness parameter k the localization border lies very low, N^. — I. The border of Breit- 
Wigner regime is = M^/48/c^ ~ 35. It means that between N^. < N < N\y Wigner 



ergodicity 



M] ought to be observed and for > N\y Shnirelman ergodicity should emerge. 



In 1974 Shnirelman lo| proved that quantum states in chaotic billiards become ergodic for 
sufficiently high level numbers. This means that for high level numbers wave functions have 
to be uniformly spread out in the billiards. Frahm and Shepelyansky 9|] showed that in the 
rough billiards the transition from the exponentially localized states to the ergodic ones is 
more complicated and can pass through an intermediate regime of Wigner ergodicity. In 
this regime the wave functions are nonergodic and compose of rare strong peaks distributed 
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over the whole energy surface. In the regime of Shnirelman ergodicity the wave functions 
should be distributed homogeneously on the energy surface. 

In this paper we focus our attention on Shnirelman ergodicity regime. 

One should mention that rough billiards and related systems are of considerable interest 
elsewhere, e.g. in the context of dynamic . 



ocalization 



localization in discontinuous quan- 



12l |. microdisc lasers Q, [l4] and ballistic electron transport in microstructures 



;um systems 

In order to investigate properties of nodal domains knowledge of wave functions (elec- 
tric field distributions inside the microwave billiard) is indispensable. To measure the wave 
functions we used a new, very effective method described in [igI]. It is based on the pertur- 
bation technique and preparation of the "trial functions". Below we will describe shortly 
this method. 

The wave functions \&Ar(r, 6) (electric field distribution E^i^r, 6) inside the cavity) can be 
determined from the form of electric field En{Rc, 9) evaluated on a half-circle of fixed radius 
Rc (see Fig. [T]). The first step in evaluation of En{Rc,9) is measurement of \Em{R^ 



The perturbation technique developed in [17| and used successfully in [17|, ll8|, ll9|, l20|] was 



implemented for this purpose. In this method a small perturber is introduced inside the 
cavity to alter its resonant frequency according to 

V - = VNiaB^ -hE%), (1) 

where v^q is the A^th resonant frequency of the unperturbed cavity, a and h are geometrical 
factors. Equation (1) shows that the formula can not be used to evaluate Elj until the term 
containing magnetic field Bn vanishes. To minimize the influence of on the frequency 
shift V — vn ^ small piece of a metallic pin (3.0 mm in length and 0.25 mm in diameter) was 
used as a perturber. The perturber was moved by the stepper motor via the Kevlar line 
hidden in the groove (0.4 mm wide, 1.0 mm deep) made in the cavity's bottom wall along the 
half-circle Rc- Using such a perturber we had no positive frequency shifts that would exceed 
the uncertainty of frequency shift measurements (15 kHz). We checked that the presence 
of the narrow groove in the bottom wall of the cavity caused only very small changes 5vn 
of the eigenfrequencies vn of the cavity \5i'n\/t-'n < 10~^. Therefore, its influence into the 
structure of the cavity's wave functions was also negligible. A big advantage of using hidden 
in the groove line was connected with the fact that the attached to the line perturber was 



4 



^ 2 



8 -(a) 
6 - 
4 





FIG. 2: Panel (a): Squared wave function |^'435(i?c) (in arbitrary units) measured on a half- 
circle with radius Rc = 17.5 cm (1^435 — 14.44 GHz). Panel (b): The "trial wave function" 
^435 (-Rc)^) (ill arbitrary units) with the correctly assigned signs, which was used in the recon- 
struction of the wave function ^435{r,6) of the billiard (see Fig. [3|). 



always vertically positioned what is crucial in the measurements of the square of electric field 
E]\f. To eliminate the variation of resonant frequency connected with the thermal expansion 
of the aluminium cavity the temperature of the cavity was stabilized with the accuracy of 
0.05 deg. 

The regime of Shnirelman ergodicity for the experimental rough billiard is defined for > 
35. Using a field perturbation technique we measured squared wave functions |\E'Ar(i?c, 
for 156 modes within the region 80 < < 435. The range of corresponding eigenfrequencies 
was from z/gg ^ 6.44 GHz to z/435 ^ 14.44 GHz. The measurements were performed at 0.36 
mm steps along a half-circle with fixed radius Rc = 17.5 cm. This step was small enough to 
reveal in details the space structure of high-lying levels. In Fig. [2] (a) we show the example 
of the squared wave function n{Rc,0)\'^ evaluated for the level number = 435. The 
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perturbation method used in our measurements allows us to extract information about the 
wave function amplitude | \l/ at ( i?c, ^) I at any given point of the cavity but it doesn't allow 
to determine the sign of "^n^Rc^^) |21|. Our results presented in [16|] suggest the following 
sign-assignment strategy: We begin with the identification of all close to zero minima of 
l"^ i\i{Rc,6)\. Then the sign "minus" maybe arbitrarily assigned to the region between the 
first and the second minimum, "plus" to the region between the second minimum and the 
third one, the next "minus" to the next region between consecutive minima and so on. 
In this way we construct our "trial wave function" j\f{Rc,6). If the assignment of the 
signs is correct we should reconstruct the wave function \E'Ar(r, 6) inside the billiard with the 
boundary condition "^^{rrjOr) = 0. 

The wave functions of a rough half-circular billiard may be expanded in terms of circular 
waves (here only odd states in expansion are considered) 

L 

"^Nir, 9) = Y^ asC^Jsik^r) sin(s^), (2) 

s=l 

where Cs = (f /o'^'""" \Js{kNr)\'^rdr)-^/^ and = 2tivn/c. 

In Eq. (2) the number of basis functions is limited to L = kNTmax = l^i""^ ■, where 
fmax = 21.4 cm is the maximum radius of the cavity. /^"^ = k^r^ax is a semiclassical 
estimate for the maximum possible angular momentum for a given kjsf. Circular waves with 
angular momentum s > L correspond to evanescent waves and can be neglected. Coefficients 
as may be extracted from the "trial wave function" '^]si{Rc,6) via 

as = ["^CsJs^kNRc)]-^ £ ^n{Rc, 0) sm{se)d9. (3) 

Since our "trial wave function" \l/jv(-Rc;^) is only defined on a half-circle of fixed radius Rc 
and is not normalized we imposed normalization of the coefficients a^: Y,s=i — ^- Now, 
the coefficients and Eq. (2) can be used to reconstruct the wave function \I'Ar(r, 6) of the 
billiard. Due to experimental uncertainties and the finite step size in the measurements of 
\'$n{Rc,0)\'^ the wave functions \I'Ar(r, 6') are not exactly zero at the boundary F. As the 
quantitative measure of the sign assignment quality we chose the integral 7 /p \'^]\f{r,9)\'^dl 
calculated along the billiard's rough boundary F, where 7 is length of F. In Fig. [2] (b) 
we show the "trial wave function" ^435 with the correctly assigned signs, which was 

used in the reconstruction of the wave function \l/435(r, 6) of the billiard (see Fig. [3]). Using 
the method of the "trial wave function" we were able to reconstruct 138 experimental wave 
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FIG. 3: The reconstructed wave function ^'435(^,0) of the chaotic half-circular microwave rough 
billiard. The amplitudes have been converted into a grey scale with white corresponding to large 
positive and black corresponding to large negative values, respectively. Dimensions of the billiard 
are given in cm. 

functions of the rough half-circular billiard with the level number N between 80 and 248 
and 18 wave functions with between 250 and 435. The wave functions were reconstructed 
on points of a square grid of side 4.3 ■ 10~^ m. The remaining wave functions from the 
range iV = 80 — 435 were not reconstructed because of the accidental near-degeneration 
of the neighboring states or due to the problems with the measurements of n{Rc,0)\^ 
along a half-circle coinciding for its significant part with one of the nodal lines of \1/Ar(r, ^). 
These problems are getting much more severe for > 250. Furthermore, the computation 
time tj. required for reconstruction of the "trial wave function" scales like tr oc 2"^"^, where 
Hz is the number of identified zeros in the measured function n{Rc,G)\- For higher N, 
the computation time tr on a standard personal computer with the processor AMD Athlon 
XP 1800+ often exceeds several hours, what significantly slows down the reconstruction 
procedure. 

Ergodicity of the billiard's wave functions can be checked by finding the structure of the 
energy surface Q]. For this reason we extracted wave function amplitudes C^^^ = {n, l\N) in 
the basis ra, / of a half-circular billiard with radius r^ax, where n = 1, 2, 3 . . . enumerates the 
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FIG. 4: Structure of the energy surface in the regime of Shnirelman ergodicity. Here we show the 
moduh of amphtudes |C^f^^| for the wave functions: (a) = 86, (b) N = 435. The wave functions 
are delocahzed in the n, / basis. Fuh hnes show the semiclassical estimation of the energy surface 
(see text). 



zeros of the Bessel functions and / = 1, 2, 3 ... is the angular quantum number. The moduli 
of amplitudes IC^'f'*! ^^'^ their projections into the energy surface for the representative 
experimental wave functions = 86 and = 435 are shown in Fig. m As expected, in the 
regime of Shnirelman ergodicity the wave functions are extended homogeneously over the 
whole energy surface [6]. The full lines on the projection planes in Fig. Hl^a) and Fig. Ht^b) 
mark the energy surface of a half-circular billiard H{n, I) = En = k% estimated from the 
semiclassical formula 



7|: ^JilT'^y -I'^-l arctan(/-y (/]^«^)2 - P) + 7r/4 = vm. The peaks 
I are spread almost perfectly along the lines marking the energy surface. 
An additional confirmation of ergodic behavior of the measured wave functions can be also 
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FIG. 5: Amplitude distribution P(^A^/^) for the eigenstates: (a) TV = 86 and (b) N = 435 con- 
structed as histograms with bin equal to 0.2. The width of the distribution P{^) was rescaled 
to unity by multiplying normalized to unity wave function by the factor A^^'^, where A de- 
notes billiard's area. Full line shows standard normalized Gaussian prediction Po{^A^/'^) = 
(l/V2^)e-*'^/2. 



sought in the form of the amplitude distribution [2^,123!]. For irregular, chaotic states 

the probability of finding the value \1/ at any point inside the billiard, without knowledge 
of the surrounding values, should be distributed as a Gaussian, ~ e~^*^. It is worth 

noting that in the above case the spatial intensity should be distributed according to Porter- 
Thomas statistics j^. The amplitude distributions P{'^A^^'^) for the wave functions = 86 
and = 435 are shown in Fig. [51 They were constructed as normalized to unity histograms 
with the bin equal to 0.2. The width of the amplitude distributions -P(^E') was rescaled to 
unity by multiplying normalized to unity wave functions by the factor A^^"^, where A denotes 
billiard's area (see formula (23) in [22]). For all measured wave functions in the regime of 
Shnirelman ergodicity there is a good agreement with the standard normalized Gaussian 
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FIG. 6: The number of nodal domains b^iv (full circles) for the chaotic half-circular microwave 
rough billiard. Full line shows a least squares fit Hjv = aiA^ + hi^/N to the experimental data (see 
text), where ai = 0.058 it 0.006, hi = 1.075 it 0.088. The prediction of the theory of Bogomolny 
and Schmit ai = 0.062. 



prediction Pq{^!A^'^) = (l/V2^)e-*'^/2_ 

The number of nodal domains K^v vs. the level number N in the chaotic microwave rough 
billiard is plotted in Fig. [6l The full line in Fig. [6]shows a least squares fit Kat = aiN+bi\/N 
of the experimental data, where ai = 0.058 ± 0.006, bi = 1.075 ± 0.088. The coefficient 
ai = 0.058 ±0.006 coincides with the prediction of the percolation model of Bogomolny and 
Schmit ^^A^/iV ^ 0.062 within the error limits. The second term in a least squares fit 
corresponds to a contribution of boundary domains, i.e. domains, which include the billiard 
boundary. Numerical calculations of Blum et al. 1| performed for the Sinai and stadium 
billiards showed that the number of boundary domains scales as the number of the boundary 
intersections, that is as V^- Our results clearly suggest that in the rough billiard, at low 
level number N , the boundary domains also significantly infiuence the scaling of the number 
of nodal domains Kat, leading to the departure from the predicted scaling Kjv ~ iV. 

The bond percolation model \2\ at the critical point pc = 1/2 allows us to apply other 
results of percolation theory to the description of nodal domains of chaotic billiards. In par- 
ticular, percolation theory predicts that the distributions of the areas s and the perimeters 
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FIG. 7: Distribution of nodal domain areas. Full line shows the prediction of percolation theory 
logio((ns/n)) = logio((s/sm,m))- A least squares fit logio((ns/n)) = a2 - r logio((s/smm)) of 
the experimental results lying within the vertical lines yields the scaling exponent r = 1.99 it 0.14 
and 02 = —0.05 it 0.04. The result of the fit is shown by the dashed line. 



/ of nodal clusters should obey the scaling behaviors: Ug oc s~'^ and ni oc l^'^' , respec- 
tively. The scaling exponents jsj are found to be r = 187/91 and r' = 15/7. In Fig. [7] 
we present in logarithmic scales nodal domain areas distribution (n^/n) vs. {s/smin) ob- 
tained for the microwave rough billiard. The distribution {us/n) was constructed as nor- 
malized to unity histogram with the bin equal to 1. The areas s of nodal domains were 
calculated by summing up the areas of the nearest neighboring grid sites having the same 
sign of the wave function. In Fig. [7] the vertical axis {ng/n) = ■^Yli^i'iT-i'^^/^^'^^ repre- 
sents the number of nodal domains n[^^ of size s divided by the total number of domains 
n^^^ averaged over Nt = 18 wave functions measured in the range 250 < N < 435. In 
these calculations we used only the highest measured wave functions in order to minimize 
the influence of boundary domains on nodal domain areas distribution. Following Bogo- 
molny and Schmit [2j], the horizontal axis is expressed in the units of the smallest possi- 
ble area s^^, {s/smin) = ^ E^Ii s / s^^l^, where sj^i = 7r(joi//c7v)^ and joi ^ 2.4048 is 
the first zero of the Bessel function Jo(joi) = 0. The full line in Fig. [7] shows the pre- 
diction of percolation theory log^o((^s /""■)) = logio((s/smj„)). In a broad range of 
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FIG. 8: Distribution of nodal domain perimeters. Full line shows the prediction of percolation 
theory logio (("-«/"-)) = "T ^°Sio((V'mm))- A least squares fit logio((ni/n) ) = as-r' logio((V^mm)) 
of the experimental results lying within the range marked by the vertical lines yields r' = 2.13ib0.23 
and 03 = 0.04 it 0.21. The result of the fit is shown by the dashed line. 



\ogiQ{{s/ Smin))i approximately from 0.2 to 1.3, which is marked by the two vertical hnes 
the experimental results follow closely the theoretical prediction. Indeed, a least squares fit 
\ogiQ{{ns/n)) = 02 — T\ogiQ{{s I Smin)) of the experimental results lying within the vertical 
lines yields the scaling exponent r = 1.99 ± 0.14 and a2 = —0.05 ± 0.04, which is in a good 
agreement with the predicted r = 187/91 ~ 2.05. The dashed line in Fig. [7] shows the 
results of the fit. In the vicinity of iogiQ{{s / Smin)) — 1 and 1.2 small excesses of large areas 
are visible. A similar situation, but for larger logi^^s / Smin) > 4, can be also observed in the 
nodal domain areas distribution presented in Fig. 5 in Ref. I2] for the random wave model. 
The exact cause of this behavior is not known but we can possible link it with the limited 
number of wave functions used for the preparation of the distribution. 

Nodal domain perimeters distribution {rii/n) vs. {l/lmin) is shown in logarithmic scales 
in Fig. [HI The distribution (ni/n) was constructed as normalized to unity histogram 
with the bin equal to 1 . The perimeters of nodal domains I were calculated by identi- 
fying the continues paths of grid sites at the domains boundaries. The averaged values 
{rii/n) and {l/lmin) are defined similarly as previously defined {ng/n) and {s/smin), e.g. 
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FIG. 9: The normalized signed area distribution NY, a foi" the chaotic half-circular microwave rough 



billiard. Full line shows predicted by the theory asymptotic limit NT,a — 0.0386, Blum et al. 



{l/lmin) = ■^Y.^Iil/l^^L where l^^l = 2-k\Js[^1/'k = 2n{joi/ is the perimeter of the 
circle with the smallest possible area s^^. The full line in Fig. [8] shows the prediction 
of percolation theory \og^Q{{ni / n)) = — y log;^o((V^mm))- Also in this case the agreement 
between the experimental results and the theory is good what is well seen in the range 
0.2 < \ogiQ{{l /Imin)) < 1-2, which is marked by the two vertical lines. A least squares fit 
\ogiQ{{ni/n)) = — t' logio((V^mm)) of the experimental results lying within the marked 
range yields r' = 2.13 ± 0.23 and = 0.04 ± 0.21. The result of the fit is shown in Fig. [8] 
by the dashed line. As we see the scaling exponent r' = 2.13 ± 0.23 is close to the exponent 
predicted by percolation theory r' = 15/7 ~ 2.14. The above results clearly demonstrate 
that percolation theory is very useful in description of the properties of wave functions of 
chaotic billiards. 

Another important characteristic of the chaotic billiard is the signed area distribution 
T^A introduced by Blum et al. [1]. The signed area distribution is defined as a variance: 
J^A = ((^+ — ^-)^)/^^; where A± is the total area where the wave function is positive 
(negative) and A is the billiard area. It is predicted [l!] that the signed area distribution 
should converge in the asymptotic limit to — 0.0386A^~^. In Fig. [9] the normalized 
signed area distribution NJ^a is shown for the microwave rough billiard. For lower states 
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80 < < 250 the points in Fig. [9]were obtained by averaging over 20 consecutive eigenstates 
while for higher states > 250 the averaging over 5 consecutive eigenstates was apphed. For 
low level numbers < 220 the normalized distribution iVS^ is much above the predicted 
asymptotic limit, however, for 220 < N < 435 it more closely approaches the asymptotic 
limit. This provides the evidence that the signed area distribution Syi can be used as a 
useful criterion of quantum chaos. A slow convergence of A^S^ at low level numbers A^ 
was also observed for the Sinai and stadium billiards In the case of the Sinai billiard 
this phenomenon was attributed to the presence of corners with sharp angles. According 
to Blum et al. jl| the effect of corners on the wave functions is mainly accentuated at low 
energies. The half-circular microwave rough billiard also possesses two sharp corners and 
they can be responsible for a similar behavior. 

In summary, we measured the wave functions of the chaotic rough microwave billiard 
up to the level number A^ = 435. Following the results of percolationlike model proposed 
by Q] we confirmed that the distributions of the areas s and the perimeters / of nodal 
domains have power behaviors Us oc s^'^ and ni oc Z""^', where scaling exponents are equal 
to r = 1.99 ± 0.14 and r' = 2.13 ± 0.23, respectively. These results are in a good agreement 
with the predictions of percolation theory ^], which predicts r = 187/91 ~ 2.05 and r' = 
15/7 ^ 2.14, respectively. We also showed that in the limit A^ — * oo a least squares fit of the 
experimental data yields the asymptotic number of nodal domains Kat/A^ — 0.058 ± 0.006 
that is close to the theoretical prediction Kat/A^ ~ 0.062 j^. Finally, we found out that 
the signed area distribution E4 approaches for high level number A^ theoretically predicted 
asymptotic limit 0.0386A^^^ [l(. 
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